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Abstract 

DNA denaturation, wetting in two dimensions, depinning of a flux line, and other problems 
map onto a phase transition with effective long range interaction. It yields giant non-universal 
critical indexes, arbitrarily large macroscopic correlation length and fluctuations at a finite distance 
from the critical temperature. In the vicinity of this region the Gibbs distribution is invalid, and 
thermodynamics must be calculated from the first principles. There are no fluctuations above the 
critical temperature. 
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Thermal unbinding (melting, coiling, denaturation) of a double-stranded DNA molecule 
is biologically important and physically unique. It yields a phase transition in a one- 
dimensional system The system is extraordinary long - the total length of a single 
mammalian DNA is 1.8m, it consists of ~ 5 billion nucleotide base pairs. Their sequence 
is related to genetic information, yet statistically it is close to a random one P]. The 
fraction of unbound base pairs as a function of temperature ("the DNA melting curve") 
is proportional to DNA light absorption at about 260nm. DNA denaturation maps onto 
a variety of other problems: the binding transition of a polymer onto another polymer, 
a membrane, or an interface j^; wetting in two dimensions 0|; depinning of a flux line 
from a columnar defect in type-II superconductors 5]; localization of a copolymer at a two- 
iuid interface 6|. DNA denaturation has been extensively studied for nearly four decades 
1 Q a fl fl HQ B a y a Q . Yet, s„.e featu..es of ..s .a..«o. we.e ove.ooUed. 
Start with its physics and model. DNA nucleotide base pairs (adenine - thymine AT, gua- 
nine - cytosine GC) are large ( "mesoscopic" ) organic molecules. Their unbinding releases 
few thousand degrees of freedom. The corresponding entropy is sks per site j3] (/cs is the 
Boltzmann constant, s ~ 10). So, while the binding (hydrogen) energy of DNA strands 
is ~ 3000°/^, DNA melts at a relatively low room temperature (~ 300°i^'), i.e. in the 
vicinity of the ground state. The Poland-Scheraga model jl| of DNA melting introduces 
the fusible AT and refractory binding energies Ei = —sksTi and E2 = —skBT2 correspond- 
ingly (Ti < T2), the boundary energy J per bound segment (J ~ 3000°/^ accounts for an 
incomplete unbinding at the boundaries), and the loop entropy — c/c^lnL per an unbound 
segment [L is the total number of nucleotide pairs there). The value of the constant c may 



vary 



mum 



from 1.5 to slightly higher than 2. Thus, at the temperature T, 
Poland-Scheraga Hamiltonian Eil^ of the adjacent bound and melted segments is related 
to the length / and the GC concentration x in the former and to the length L in the lat- 
ter. Calculated from the energy —sksT per site of a completely melted DNA (T is the 
temperature), 

EiLx = skslST + J + c/cfiTlnZ/ — skslix — x)AT. 
6T = T-f, f = Tix + T2{l-x), AT = T2-Ti, 

where x is the AT concentration at an entire DNA. Parameters T ~ 310K, AT ~ 407^^ de- 
pend o. tKe DNA .o.ution Q. TKe P„.a„d-ScKe.aga n„de. © allows fo. a .t.a,-Mfo..wa.d 
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calculation of the homopolymer thermodynamics. The calculation yields macroscopic fluc- 
tuations. I prove that heterogeneity just renormalizes the interaction parameters c and J . 
The renormalization depends on temperature and heterogeneity. This implies an unusual 
non-universal transition, which invalidates the Gibbs distribution in its close, yet macro- 
scopic vicinity. Start with the well known case of a homopolymer 0, Q, 0- There x = x, 
the last term in Eq. is missing, and Eil^ = E{l,L) depends on / and L only. Then an 
entire Hamiltonian H = E{1^"\ L^"^^) describes an ideal gas of the pairs It 
relates the free energy / per site to the normalization condition for the Gibbs probability 
piL of given / and L: 

oo 

piL = [exp-{l + L)<f)~E{l,L)/kBn 5^P/L = 1, <p=-f/kBT (2) 

l,L=l 

When < expi-J/keT), Eqs. (PH) yield 

/oo 
exp{-L(f))L-^~^'dL = (0 + r) exp{J/kBf); r = s5T/f; ci = c - 1 ^) 

Consistent with the Landau- Peierls theorem for the Hamiltonian when Ci > 1, Eq. ()3b|) 
yields phase transition. Then, by Eq. Q, = does not allow for any excitations of 
a completely melted polymer. This is speciflc for the Hamiltonian which depends on InL 
only — when L = oo, any excitation implies an inflnite energy increase. Dependence on 
In L yields other unusual implications also. Transition is non-universal - its critical indexes 
depend on ci. Immediately below the critical temperature 0,0,Q,S| Tc, 

= ^ifci>l; 0~ -^In^ if ci = 1; ~ ^^/^Mf 1 > Ci > 

(3) 

e = {n- T)/{T, - f ); r, = (T, - f )/f = (1/sci) exp(- J/fc^T); 

As anticipated, the critical ci = 0, while J/ksT ~ T/ AT ~ s ~ 10 implies, by Eq. Q, a 
very narrow width of the transition ~ [T^ -T)/T ^ 10"^ (i.e. - T ^ IQ-^K ), and its 
very close proximity to the ground state melting temperature T. Once the free energy (jSJ is 
known, the Gibbs probability (0) allows one to calculate any thermodynamic averages and 



fluctuations. The average (denoted by a bar) relative number uj = l/L oi the bounded sites, 
which is measured via light absorption, the average length L of a melted segment, and their 
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relative mean squared fluctuations Auj/u, AL/L are: 

o) ~ ci exp{J/kBT) ^ 1; Au/uj ~ 1 
I ~ 1 if ci > 1; I ~ 0'^"^ if ci < 1 

AL/L ~ 1 if ci > 2; AL/Z ~ /-Sci-i if 2 > ci > 1; AL/L ~ if 1 > ci > 

(4a) 

Thus, Auj/u, AL/L are never small, while AL/L oo when T ^ Tc and ci < 2. A more 
physically meaningful fluctuation is 



A*uj/oo=\uj-io\/uj r^uj-^^ <^1; A* L/ L = \L - L\/ L ^ I (4b) 

It demonstrates, in particular, that a characteristic \L — L\ ~ L implies a characteristic 
I InL — lnZ| ~ 1 i.e., <^ InZ when Ci < 1 and L oo. 

Consider heterogeneous DNA. When temperature increases from T to T+(5T, the Poland- 
Scheraga Hamiltonian complements the energy increase of an "average" bounded segment 
(the first three terms) with the energy decrease of a refractory bounded segment (the last 
term). I prove that in the vicinity of the DNA melting temperature, the last term may 
be replaced with its thermodynamic average for given lengths of the successive bound and 
unbound segments. (Such replacement is equivalent to an unusual mean field approxima- 
tion, which becomes accurate at the phase transition and which technically reduces to a 
constrained summation in the partition function). The resulting Hamiltonian describes a 
homopolymer with the renormalized loop entropy. The renormalization, and thus the phase 
transition singularity it determines, are non-universal and depend on the DNA parameters. 

Physics of DNA melting was elucidated in Ref. 0. (All following statements are later 
accurately verified). The segments, rich in the fusible AT, melt first, while the richest 
in the refractory GC melt last. When c > 1 in Eq. (0), bounded segments completely 
vanish at a finite critical temperature Tc- There L ^ oo and the effective boundary energy 
(J + cksT In L) cxD if c < 2. Then, the excitation energy also — > oo, DNA approaches its 
ground state, fluctuations of I and x vanish, and the length of a ground state (i.e. sufficiently 
refractory) bounded segment is p] oc InL ^ oo, to compensate the effective boundary energy 
in Eq. (PJ). Sufficiently close to Tc, / exceeds any finite correlation length, and the probability 
w{l,x) of a given x at such / is Gaussian. Since fluctuations of / and x vanish at T^, it is ~ 
the thermodynamic probability //L of a bounded site. So, I and x yield L{l,x) ~ l/w{l,x). 



Thus, physics of DNA melting suggests, and further calculation verifies, that the values of 
/ and L ^ / at a given temperature determine the corresponding value of x according to 
w{l, x) ~ 1/ L. The Gaussian w implies v^(x — x) oc \J\n\i/w{l^ x)] ^ 1 where the factor in 
w ^ l/L may be disregarded with negligible error. Thus, x — x in Eq. (0) may be replaced 
with its thermodynamic average according to 

l/L = {l/27iD^y^^exp{-u^y, = l{x-xf/2D'^- D"^ = x{l - x) j^) 

In fact, large J/ksT ~ T/AT ~ s ~ 10 allow for Eq. (jKEj) already slightly above T. 
Equation (^, complemented with the unusual mean field approximation ()5bp for x, yield 
the renormalized Hamiltonian E*{1, L), which depends on the variables I and L only. In the 
leading (in l/L <^ 1) approximation it equals 

E*{1, L) = slkeST + J + ckeT In L - skBDATV2l InL. (5) 

The last refractory term accounts for the thermodynamic average of x for given values of /, L, 
X and AT in the Poland-Scheraga model for a heteropolymer. The Hamiltonian describes 
a "renormalized" homopolymer, and Eq. (j2I), where E{l,L) is replaced with E*{l,L), yields 
its exact free energy. The knowledge of the free energy allows one to calculate all averages 
and their fluctuations, to estimate the accuracy of the crucial approximation ()5b|) . to prove 
that its inaccuracy vanishes at the phase transition, and thus Eq. (0) accurately determines 
transition singularities. It also presents an exactly solvable model of a heteropolymer and 
DNA. The competition in Eq. ^ of the energy increase and decrease, correspondingly in 
the "average" first and last "refractory" terms, yields a high and relatively narrow E*{l,L) 
minimum at the ground state I = 1^ = 0.5{DAT/6T)'^\nL (which is indeed oc InL as stated 
earlier). The expansion of E*{l,L) in (/ — Im) non-universally decreases the factor c in the 
loop entropy by s{DAT)^ /2T6T, and Eq. (j21), with E replaced with the expanded E*, after 
a straightforward calculation, yields 

/oo 
{lnLy/^L-^-^exp{-(f)L)dL = M (6) 

where 

5 = ci - 7; 7 = s{DAT f/2T5T- 
M = 7i-^/\2ci)-^/\sDAT/ffexp{J/kBf) > 1. 
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Note that the left hand side of Eq. © depends on and S only. Thus, Eq. © reduces five 
dimensionless parameters {J/T, AT/T, T/T, x, c), which determine in a non-renormalized 
case, to two parameters (5 and M). When ci < 1 and 4> <^ Tc — r, Eq. ^ maps onto Eq. 
(j3bp . where Ci is renormalized into the temperature dependent 6 (which, unlike Ci, may be 
of any sign and which dominates the temperature de pen dence in the vicinity of the phase 



transition). By Eq. > 0. Since ci > 
at finite temperature T = Tc. There 



ature ae pen aenc 

!ilQ,Hliil,Q, 



is achieved (as stated) 



5Tjf ~ (s/2ci)(DAT/f )2, 5T, = T,-f 



(8) 



Note that, by Eq. ©, ~ 3°K. When - T <^ 6T^, then 6 = ci{e* - 6), where 
6 = {Tc — T)/6Tc is the relative distance to the critical temperature. By Eq. ()4b|) . 5 -C 1 
implies 6lnL ~ 6lnL, and thus verifies the derivation of Eq. (jSj). At Tc, by Eq. (jHJ, 
6c ~ 0.01, i.e. it is very close to the critical Ci = — c/ Eq. Q. By Eq. L oc l/cp oo 
when T — Tc. This, and / ~ oc InL verify all previous estimates. When 5, ^ 1, 
asymptotics in Eq. (jHl), where M 1, yield an unusual non- universal singularity: 



~ [(30/47r^/2^*)v/ln(^V^)]'/"i^*; when ^ < r (9a) 



~ [(2^77r0)^ln(0/0*)]^/2cie. ^^^^^ ^* < ^ < 1 (9b) 

Consider the implications of Eqs. (f7ti9b|). In natural DNA J/ksT ~ T/AT ~ 10, 
D ~ 1/2, Ci ~ 1. So, in the immediate vicinity of Tc, where 9 < 9* 0.01, the order 
of the transition, by Eq. (j9a|) . is l/ci9* ~ 100, i.e. giant. The order is non-universal, it 
depends on the DNA parameters Ti, T2, x. The values of Ti, T2 depend on the hgands and 
their concentrations in the DNA solutions 0,0, which may be manipulated experimentally. 



Non-universality in Eqs. ()9al I9b|) is related to the competition of the refractory and loop 
entropy terms in Eq. (jSj), which renormalizes the loop entropy, and thus the singularity. 
The width of the transition is very small, yet macroscopic. The crossover from Eq. ()9b|) 
to Eq. (jnil) occurs when (Tc - T)/Tc ~ 10"^ Then Tc - T ~ O.OlisT (cf 5Tc = Tc - 
T ~ 3K). In the approximation of Eq. ©, the probability density Pl of a given L is 
= M-i(lnL)i/2L-i-^exp(-0L). So, by Eqs. (0 Z oc 1/0 oc exp[l/(Tc - T)] 
exponentially increases to L ~ 10^'^ at the crossover. Thus, even in a solution with ~ 10^^ 
DNA nucleotide base pairs, all DNA molecules completely unbind in the interval (j9bjl . So, 



at a small, yet macroscopic distance ~ O.Olii' from Tc, the effective long range interaction 
exceeds any macroscopic size of the system. The system can no more be divided into weakly 
interacting subsystems, thus the Gibbs distribution is invalid, and thermodynamics must be 
calculated from the first principles. The fraction of bounded sites is correspondingly small 
there, and the observably quantity is the temperature of complete melting of a finite DNA. 
If its length is A^, then L = at the temperature T/v, where 

^^ = (T,-T^)/Te~l/lnAr. (10) 

The mean fiuctuation A9n- of may be estimated from L{9n + A9n) — L{9j^) = A*L{9j\f). 
Similar to Eq. (gEI), A*L ~ L and thus Q 

A9n/9n l/lnN. (11) 

Such fiuctuation is macroscopic and easily observable. In DNA this situation is related to 
mesoscopic size of base pairs (which yields large J/ksT ~ s ~ 10, thus small 9*), and to 
DNA heterogeneity. By Eq. heterogeneity effectively replaces fixed ci with 6. The 
latter decreases to 6c ^ I (at T = Tc) and scales the transition order with l/6c in Eq. 
and with 1/5 in Eq. ()9b|) . This may be characteristic of any sufficiently strong long range 
interaction. 

By Eq. ()10|) . natural DNA always yields 9 9* , i.e. the essential singularity ()9b|) in 
9 (X Tc — T. (This was predicted in ref. ^). By Eq. (jH)), it proceeds in the interval 
Tc — Tr^ O.OlTc ~ "iK. Sufficiently close to T, the length L may reach the correlation length 
of the sequence. Then the distribution w{l,x) becomes non-Gaussian. This alters Eq. (jSJ 
and the melting curve 4>{T). 

Below T DNA is mostly bounded, and only anomalously fusible segments melt. Their 
probability yields the equation which replaces Eq. (jSj). Their melting proceeds in an entire 
interval AT. Until sufficiently high temperatures, when the number of segments, which 
melt nearly simultaneously, becomes large, the DNA melting curve exhibits their successive 
melting. It is explicitly seen in experiments 0,01 • Thus, in a general case there are three 
distinctly different temperature intervals: 9* ~ 0.01, i.e. Tc — T 0.03K; 6' ~ 1, i.e. 
Tc~T 3K] and AT ~ 40K. 

A giant order transition (jUaf) may be observed only when the total number A^ of base pairs 
is much larger than L at the crossover to Eq. ()9b|) . This implies In A^ > 1/9*. Since, by Eq. 



(jHl), 0* ~ 0.005^-^/3^ gQ D j^^g^ ^ 0.03(lnA^)3/4. Qn the other hand, the derivation of 
Eq. (jni) imphed the large renormahzed term. At the crossover this means D > 0.03. In the 
interval 0.03 < D < 0.03(ln A^)^/^ non-universality of the giant critical index in Eq. ()9ap may 
be studied (e.g., via its dependence on AT, which changes together with the concentration 
of solvents in DNA solution 7]). 

Presented theory may be numerically tested. Once the ground state is accurately deter- 
mined analytically ^, computer simulations allow for the study of its fluctuations. 

The approach is applicable to other problems also. 

To summarize. DNA unbinding with temperature proceeds from piecewise melting of 
fusible domains, to essential singularity, to giant (~ 1/6* > 100) order phase transition. The 
latter may be observed when the AT or GC concentration is between 0.03 and 0.03(ln A^)'^/^, 
where A^ is the total number of nucleotide pairs. In the vicinity of complete melting the Gibbs 
distribution is invalid, and thermodynamics must be calculated from the first principles. 
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